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On a Form of the Solution of Laplace's Equation Suitable for 

Problems Relating to Two Spheres, 

By G-. B. Jeffery, University College, London. 

(Communicated by Prof. L. N. G. Filon, F.R.S. Eeceived April 11, and in 

revised form June 10, — Eead June 13, 1912.) 

§ 1. Introduction. 

The problems presented by the motion of two solid spheres in a perfect 
fluid have been attacked by various writers * In each case the method has 
been that of approximation by successive images, and it appears that no 
general analytical method of solution has been developed as in the case of 
the analogous problems for the sphere, ellipsoid and anchor-ring. In 
this paper a general solution of Laplace's equation is obtained in a form 
suitable for problems in which the boundary conditions are given over any 
two spherical surfaces. A similar solution is obtained of the differential 
equation of Stokes' current function. With the aid of these results it is 
theoretically possible to determine completely a potential function when its 
value is specified over any two spheres. The method is illustrated by an 
application to the electrostatic field of two charged conducting spheres. In 
this case the method leads to a simple expression for the capacity of either of 
the spheres. 

The co-ordinates employed are defined by rotating about the z axis the 
system of circles, in any plane, through two fixed points on the axis and the 
orthogonal system of circles. Thus, if x, y, z, are the Cartesian co-ordinates and 
p = ^/(x 2 + y 2 ), and the distance between the fixed points is 2a, we have a 
system of orthogonal curvilinear co-ordinates u, v, w, where 

• i b ~\~ % (z -f~ ct) ti 

u+iv = log ^ .) — {, ^ = tan~ 1 ^. (1) 

3 p + i(z — a) x v ' 

The surfaces u = constant will then be a series of coaxal spheres having 
a common radical plane u = 0. It is obvious that the origin and the value 
of a can be so chosen that any two given spheres will be included in the 
system. These co-ordinates are similar to those employed by Hicks in his 
memoir on " Toroidal Functions/'f the difference being that in the present 

* R. A. Herman, " On the Motion of Two Spheres in a Fluid," ' Quart. Journ. of Maths.,' 
1887, vol. 22 ; W. M. Hicks, "On the Motion of Two Spheres in a Fluid," 'Phil. Trans.,' 
1880 ; A. B. Basset, " On the Motion of Two Spheres in a Liquid," ' Proc. Lond. Math. 
Soc.,' 1887. 

t 'Phil. Trans.,' 1881, p. 609. 
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case the circles are rotated about the line through the limiting points instead 
of about their common radical axis. Further, in one case the surface 
conditions are given over spheres, while in the other they are given over 
" tores " or anchor-rings. For these reasons, although that part of Hicks' 
paper which deals with the general theory of curvilinear co-ordinates has 
been of very considerable help to the present writer, the functions obtained 
and the methods suitable for their applications are very different in the two 
cases. The following results w T ill be of use ; they are easily obtained from (1) 
and are set down here for reference : — 



a&mv a sinh u /ox 

cosh u— cos v cosh u — cos v 



'du\ 



2 



If Ss U) Ss v denote the lengths of elements of the curves u = constant and 
v = constant respectively, 



Bs u = (M + MVbv = -^ Bv. 

C\uuJ xovl J cosh u— cos v 

Bsv = ( (h) 2 + (ffy Bu = — ^ Bu. 

LVow \ov/ J cosh u — cos v 



> (4) 



If r be the radius of any sphere of the system, and d the distance of its 
centre from the origin, 

r == a J | sinh u | , d = a coth u. (5) 

§ 2. Solution of v 2 <£ = 0. 
It is well known that if u, v, w, be any system of orthogonal co-ordinates 
Laplace's equation may be written in the form 



V 2< £ = h^h 



a (h^ d$\ a_ / h 8j>\ + ±/fo fy 



hu \Ji 2 hz du j dv \hjii dv j dw \hih 2 div)_ 



0, (6) 



In the present case, ^, v, being conjugate functions of p, z. 

'"•=|J + (I'=(|WI)'=^ 

and if w = tan" 1 2, hi = -=. 

x p* 
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Substituting these values in (6), we have, as the form taken by Laplace's 
equation in u, r o i w co-ordinates, 

In order to solve this equation, write <£ = tyl\/p. We obtain 

du>* do> + 4 P 2 IW U/ J P L\9p/ + W J 9w 2 ~ ' 
which hecomes, in virtue of equations (2) and (3), 



8^ 2 * 3v 2 sin 2 ^14 dw 2 / 



(8) 



After the usual manner we will seek a solution of the type -\/r = U, V, W, 
where U, V, W, are functions of u, v, w, respectively. Substituting in (8) 

we at once obtain 

.= cos (mw -fa) 

1 d 2 U 1 d 2 V , l-4m 2 , , , ivo 

Hence U = A cosh (% + -J) u + B sinh (w + 1) w, 

and ^ +((^ + i) 2 -^^ : ) V = 0. 

cfe 2 \ 7 4sm 2 W 

Put V = >^/(sin v) x an( i we obtain 

^-\-cotv-%' + (n 2 + n — ^-- )v = 0. 
dv 2, dv \ sm J vl 

Finally, writing cos v = ^ we have 

the solution of which is well known to be 

where P„ m 0*) = ( 1 - n 2 )™' 2 ^J"^ , 



P n (p,) being the Legendre function of order ». Q B m (jj,) is the corresponding 
function of the " second kind." Hence 

V = v^sin v) {aP„ m (>) + &Q» m (/*)}, 
and a particular solution of (7) may be expressed in the form 
(j> = -(/'(cosh u—fi) cos (mw + a) {A cdsh (n + ^) ti + B sinh (n + -|) %} 

{aP» m 0*) + &Q»»0*)}. 
We shall find that for physical applications it is sufficient to confine our 
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attention to integral values of n. Moreover, the solution corresponding to 
?i = — (v+1) is identical in form with that corresponding to n = +v. It 
will, therefore, be sufficient to consider only positive integral values of n, and 
since by the nature of the function P w m (/a), m^ \n\, we may write the solution 
of (7) in the form 

m—ao n — co 

cf> = ^'(cosh^ — yLt)2 Xcos(mw + ct) {A min cosh(n + ^)ii + B mtn smh(n + ^)u} 

§ 3. The Properties of the Function </>. 

The function <j> and its first differential coefficients must be finite and 
continuous at all points of the field except those which correspond to some 
special physical condition such as a source or charge. P n m (cos?;) is finite and 
continuous for all real values of v, but Q n m (cos v) becomes infinite when 
v = or 7r, and hence cannot occur in an expression for <f> which is to hold 
throughout a region including any part of the axis of z. It is to such cases 
that we confine our attention in the present paper. At the limiting points of 
the system u becomes infinite, and hence the forms of </> appropriate to regions 

including the points u = _ are respectively 

to = qo 7i=co 

<£ = ^(coshw — fi) 2 t A^ w cos(mw + a)e T < n+ *> M P n TO (/4 

It will be convenient for reference to have the general value of dcf>/du. 
Differentiating (9) with respect to u we have, omitting the Q terms, 

Jr = — — ^ ^ cos / ryw _|_ a \ (£ sinh % cosh (n -f A) % 

dw 2 v /(cosn%— /A) m=0 »=m 

+ B m , w sinh w sinh (?i + -J) u + (2 w + 1) A m> w cosh % sinh (% + i) u 
+ (2 7i + 1) B TOi TC cosh ^£ cosh (% + A) u 

— .(2 w + 1)/a [A WtW sinh (^ + J) w-f B m§ w cosh (?i + J) ^]}P« m (/a). 
Applying the recurrence relation 

(2n+l)fiV n m (fji) = (w + m)P?_iOA) + (n-m + l)P»+iO*), 
we have, after further reduction, 

§£ = i v 

-m.=0 n=m 

x sinh (w + •§•) ^ + [^A m> n —(n—m) A m> w _ x ] sinh (n — J) % 

+ [(ti + 1) B TOj w — (ti + m + 1) B m? n+1 ] cosh (ti + f ) % 

+ [>B min — (ra— m)B TO , n _i]cosh(7i— %)u} P„™(/a), (10) 

provided that (9) may be differentiated term by term and rearranged. 



m = oo n=<x> 



*% = J (cosh w— fi) * 2 S cos(m^ + a){[(7i + l)A, w>n — (7i+m + l)A TO|W+1 ] 

OW m=0 n=m 
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The coefficients in the function <£ are not in all cases identical in the two 
regions into which the field is divided by the plane u = 0. Let (9) be the 
expression for <£ in that part of the field for which u > while a similar 
expression with accented coefficients holds throughout the other half of the 
field. In order that <£ may be continuous across the plane u — 

A' —A 

while, in order that d<f>/du may be continuous, we must have (provided that 
d^/du is given by (10)) 

(n— l)B TOf „— (n + m + l)B Blfn+1 +%B in|tt — (n— m) B OT|n _i 

= {n + l)W m , n — (w-f-m + l)B / m)n+1 + ^B / mjn — (?i— m)B / mtW _ 1 

where w = m and B m , w = if n < m. 

(w + m + l)(B Wift -B IB , B+ ])-(?i-m)(B M)M . 1 -B W|n ) 

= (w + m+1) (B^ — B' m , n+1 ) — (w— m) (B' m>n _i — B' m , n ). 

(?l + m+l)(/3 mjW + 1 -/3' m ^ + l) = fa — m)(£ TO| „ — £',*,„). 

This relation holds for all values of n such that n = m. Hence 

Pw-fi'm t n = if w = m + 1, 
which at once leads to the result 

J? _ p ' ~R — TV 

UnijU iy m,n — -"^ra 1J m,m* 

These forms occur whenever the field includes doth limiting points together 
with any part of the plane u — 0, for any expression for $ of the form (9) 
which is finite at one of these points becomes infinite at the other. 

Eor a unit charge at the points u = _ °° the potential is given by 

<£ = { p 2 + (z + a) 2 }-% = ~k — s/ (tosh. u-fi)e ±u l 2 . 

Thus, as in the case of spherical harmonics, the harmonic of zero order 
corresponds to a point charge. The analogy, however, extends no further ; 
the nt\i harmonic does not correspond to a charge of the wth order. It may 
be shown that, if n be the radius vector from the poie u = + oc , 

nl gpr(±) = v/(cosh («-,*) [«*---C ie i»p 1 0*) + -O 8e s»P,( At ) 

+ ... + (-l)™C/ + «»P r ( /i ) + ...] 
where "G r is the ordinary binomial coefficient. 
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§ 4. Case of Symmetry about an Axis ; Stokes' Current Function. 

When the field is symmetrical about the axis of z, we obtain the form of 
the function <j> by putting m = in (9) 

■71 = 00 

=r ^/(cosh^— fi) % {A„cosh(^ + J)w + B»sinh(w + J)^}Pn(^). (11) 

Many physical problems are more easily interpreted by means of " lines of 
force," or " lines of flow/' than by the potential. In the symmetrical case 
there exists a function tJt corresponding to Stokes' current function, such that 
the surfaces yfr = constant are orthogonal to the equipotential surfaces, and 
which is connected with <£ by the following relations* 

dp dz ' 3# 3/5 * 

Transforming to u, v co-ordinates, 

f dyfr du dyjr dv\ ___ /^^i^^ 

9w dp dv dp) \du dz dv dz, 

f dty du dty dv\ __ _ /3j> die d(j> dv\ 
x du dz ~ch dz) \du dp dv dpi 

t» n i ,„..,,,. 3% 3v -, 3?^ 3v 

By a well-known property of conjugate functions ^- = ^- and =- = — *-. 

Hence 

du fd^r __ d<f>\ diold^r 3<£ \ __ q 
3p \3w- cW 3s; \3# 3^/ 

^^-0^ + — f^ + P^ = 0. 

02 \3w <W dp\dv du/ 

Hence, since ( ~-l ■+[p-j r 6 u > 

H - Z±. d±-_ d£ (12) 

3£6 3v dv. ow 

When a potential exists, and is continuous, we may write 

iL/s^—.?. /?£\ — o (i3) 

du\dn/ dv\du) 
Substituting from (12), we obtain the differential equation satisfied by ^, 

3iAp du) dv\p dv ) 
an equation which is somewhat similar in form to Laplace's equation and 

* Lamb, ' Hydrodynamics,' p. 117. 



\auj \dvl _ 



= 0, 
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which yields to an analogous mode of treatment. Writing y}r = ^/p.H, we 
have 

Su 2 + dv 2 2p \du 2 dv 2 ) 4p 2 
which becomes in the case under consideration 

^n + 85T 3n_ = 

du 2 * dv 2 4sin 2/ y 
Assuming, as before, a solution of the form II = UV, 

1 d 2 JJ 1 d 2 V , 3 , , / , i x2 

Tf 71 = — tt Ti+J-^T" = con stant = O + i) . 

Hence U = a cosh (% + J) w + 6 sinh (n + J) ^ 

In order to solve for V, substitute V = —% — . We obtain 

v /sm / y 

-7-4 — cot v-z — \-n(n+l) v = 0. 

Now write fi = cos v and the equation becomes 

(W)|*+»(*+l)% = 0, 

the general solution of which is 

X = A{P„-i(/a)— P n+ i(^)} + B{Q»-i(/A) — Q»+i(/i)}. 

As for the function <£, it is only necessary to consider positive integral 
values of n, and the solution of (13) which is finite on the axis is 

1 ft — CO 

^ = — n? {a n cosh (% + £•) & + & n sinh (% + AW} 

v/(coshw — /*)*=<> j 

{Pn-l(/^)-Pn + l(/^)}. (U) 

§ 5. Properties of the Function yfr. 
It is obvious that the remarks made in § 3 on the continuity of the 
function <j> apply equally to the function -^. Differentiating (14) with respect 
to u and applying the recurrence relation 

a* [P.-iG*)-P, + iG0:i = ^zi [ p »-»(A*)-P»(A*)] + ^|[P» W-P.+. G*)l 

we have 

*r = 777 — u vr 2 -77 — rV 1 — 7-?[( 2 ^ + 1 )^n— (2^+3)^ tt+ i] 

d% 2 (cosh u—fi)? n=z0 2n + l [_n + l K v 7 J 

x sinh (% + f) ^ 
— [(2%— l)^ n _i— (2w + l)a n ] sinh(?i — £)w 

+ — - T [(2^+l)&„— (2ti + 3) & w+ i] cosh (» + |)w 

%+ 1 

-i[(2»-l)J w _ 1 -(2» + l)5 )l ]oosh(»-i)«j. [P»-i(/*)-P,.+i(m)]. (15) 
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As for cj> the coefficients a must be the same on both sides of the plane 
u = 0, while to secure the continuity of the potential gradient we must have 
provided that (14) may be differentiated term by term 

(Ai + l — fin) (fin — Bn-l) = r~r (fin+1 — ' fin ) (fin — fin-l)> 



n+1 n n+1 n 

where fi n = (2n+l)i n . 

Putting n = 1, 2, 3, ..., n—1, and adding, we have 

~(fin-fin-l-fin+fin-l) = fil-fio+fil~fio' = 2<>, 

where c is a constant. 

Multiplying by w and again putting w = 1, 2, ..., n, we obtain on adding 

n 

n -0 n ' = o -0 o ' + 2c^i, 

i 

A simple example of this form is afforded by the case of a uniform field 
parallel to the axis of z in which 

, ,., __ 2n(n + l) 

u n u n — o ~z — • 

2n + l 

Relations between the Coefficients in <f> and ty. 
Let 

<£ = (cosh^— jul)%% {A n cosh(n + J)^ + B w sinh(yi + J)^} P w (//,), 

yjr = (cosh^ — //,)"^ {»» COSh (ft, + J) U + 6 W sinll (tt + 1) %} {P w -l(/i) — P w +l(yL6)} ? 

be functions specifying the same field. 

The coefficients will be related in such a way as to satisfy equations (12), 
which for our present purpose may be written in the form 

*hk. = ^(/^ 2 — 1) 8<fr . 8i/r __ a 3$ , 1 fi . 

3% COSh% — /i,3//, ' 3yL6 COSh U — fJb du ' 

Differentiating the expression for <f> with respect to /x and applying the 
well-known recurrence relations we obtain 

a \ ~~ I J? = — |a (cosh w— //,)"*£ ^ , ' {(A w — A n+ i)cosh(n+|-)^ 
cosh^— //, d/x L 2^+1 ^ ^ 

+ (A»— A»- 1 )cosh(w— J)?*-f (B a — B B+ i)sinh(w + |-)^ 

+ (B„— B n -i) sinh(^— %)u} {P w ~i(/*) — Pn+i (/a)} 

Identifying this expression with the value of dyjrfdu given by (15), we 

have 

(2n + l)a n —(2n — l)a n -i = — an(B n — B„-i),^ 

r C 1 ?) 

(2^ + 1) b n —(2n—l)b n -i = — a7&(A» — A»-i). J 
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It is easily verified that these are also the relations necessary and sufficient 
to satisfy the second of equations (12). Putting n = 1, 2, 3, ..., n, and 
adding, we obtain 



n— 1 

(2n+l)a n = — ctnB n + a % B m 

o 

(2?i + l) b n = •— <mA w -f a £ A 



>■ 



o ^ 



(17') 



§ 6. 7%e Potential and Capacity of two Spherical Conductors, 

We will conclude the present paper by an application of this method to a 

classical electrostatic problem. Let it h u 2 be any two spheres such that 

u\ > but u 2 is unrestricted. The potential <£ is constant over each of these 

spheres, and we can without loss of generality suppose it to be zero over the 

surface u 2 and V over the surface tti. It is obvious that $ will be of the 

form 

(f> = (cosh% — ^)*SA n sinh(?i + i)(^— u 2 )V n (fi) 9 

from which we have 

V (cosh u—fjb)~-% = XA n sinh (n+%) (ui—u 2 ) P n (/*). 

The left-hand side may be written 



71=00 



v /2.V^-^(l-2/>t6- M i + ^- 2M i)-i = ^/2. VSe" (B+i) ^P„(/*), 

71 = 

since u\ > 0. 

Equating the coefficients of P«(aO, we obtain the value of A m and the 
potential function becomes 

r v v *=oo smh(w + i)(wi— ^2) v 

If cri, cr 2 denote the surface densities of the charges on the spheres u h u 2) 
respectively, we have, by Coulomb's theorem, 

a — _3^ _ 8^ _ cosh^i— cos'y 3^ 

4fc 77*0*1 — ^ — #*> ' ;=: 

on os v a on 



u=u. 



= -¥_ (cosh } te - {n+ $ Ul r sinh ^ — 

y/2-a n=0 Lv / ( COSn ^l--/ u ') 

+ (2n+l) -/(cosh % — ju,)coth (w + J) (^1 — ^2) ^Pn(Afc), 
4. — cosh 1^2 — cos ^ 3$ 



71=00 



■(cosh ^-j*)* 2 ^Tirz: / iw<t , — tt p »0)- 



\/2.a ^ n==0 sinh (n +|)(% — ^2) 

Using Maxwell's notation we will denote the coefficients of capacity by 
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%\\, $22 and the coefficient of induction by qi2. To obtain these we put V = 1 
and integrate the surface density over the appropriate sphere. The following 
integrals will be required, 

_ 1 (cosh u— jj^p r 2 Ti + 1 



and 



'+i P (a) e -(n+i)\u\ 

_ 1 (coshw— ya/ ^ v smh |% 



^ii = 27T fpcricfe 



w 



+ln =°° / M+ iw, f sinh^i 



+ (2^ + l) °°^< n + *)^7 ^> P. 00^ 

x ' (cosh%— /x/ J 



71 = 00 



= a t e~( 2n+ V u i [1 + coth (n + %) (ui—u 2 )] 



n = oo 



2a 2 {e^n+ljiti^^n + l)^-^ 

w=0 



2i'2 = q.21 = 2nr\p<T2cls u 



,^.j 71 = CO 

2 



(2^ + l)<r (n+i)M > 



2^/2 . a J _ x w= o (cosh %2 — //r)a sinh (w + 1) (%i — 1*2) 



V n (/a) dp 



71=00 



= -^2 — 

ws=0 sinh (n + £) (^1 — %i 2 ) ' 

71=00 



if ^ 2 < 0. 



71=0 



From considerations of symmetry it is clear that we should obtain for the 
third coefficient 



7l = QO 



q 22 = 2aS{r^ + - 1 )^-r^ +1 )«<}- 1 . 



* 



n=0 



* These expressions are not those obtained by Maxwell (' Electricity and Magnetism 
3rd ed., vol. 1, p. 271), but it is possible to transform them into the forms there 
given, e.g. — 



n 



_oo 



q n = 2a 2 [ei^+i)^- e (2n+i)uq-i 

7i=0 

n = 00 r = 00 
= 2<22 2 e~(r + l)(2n + l)u l +r(2n + l)u 2 , 

7i = r=0 

Since w 2 <0<^ 1 the double series is absolutely convergent, and we may interchange the 
order of summation 

= 2^2 2 e-( 2 r + l)(n + l)ih + (2r+l)nu 2 
n = co 

= a2 cosech \n {u x — u 2 ) + u x 1, 

n=zQ 

which is Maxwell's form. 
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When U2 > one sphere is enclosed within the other and there are slight 
differences in the coefficients due to the different forms taken by the integrals 
in these circumstances. We obtain 



71=00 



2ii = -212= -22i = 2aS{«( J, » +1 ) M i-^ +1 )«»}- 1 , 



71 = 00 



222= •? + 2aS{g<g» +1 >*-gP»+ 1 >*}- 1 . 

sinh% «=o 



For two equal spheres % 2 = — -%. 



W=oo 



2n = 222 = ^2 cosech (2^+1)^1; 

w=0 



W = oo 



q l2 = — a Se"" (2n+1)M| cosech(2^+l)%. 

w=0 

In the case of a sphere in the presence of an infinite plane, u 2 = 0. 

j n = -2i 2 a= 2&2{^ +1 > w i--l}~ 1 . 
We have tabulated below some of the more important of these functions. 



Table 


I. — Capacity of Two Equal Spheres. 


u 


2r __ diameter of spheres 


?H. 


w • 


d centre distance 


r 


0-2 


-9804 


1 -8368 


0-4 


-9268 


1 4796 


0-6 


'8442 


1 -3055 


0-8 


'7526 


1 -2011 


l'O 


-6508 


1 -1284 


1*2 


'5526 


1 -0916 


1-4 


-4654 


1 -0609 


1*6 


-3888 


1 '0408 


1-8 


'3224 


1 '0273 


2-0 


'2662 


1 0183 


2-5 


-1630 


1 -0067 


3-0 


'0996 


1 -0025 
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Table II. — Capacity of a Sphere in the Presence of an Infinite Conducting 

Plane. 



u x . 


r _ radius of sphere 


2'n > 
r 


p distance of centre from plane" 


0-2 
0*4 
0-6 
0-8 
1-0 
1-2 
1-4 
1-6 
1-8 
2-0 
2-5 
3-0 


-9804 
0-9268 
-8442 
0*7526 
-6508 
-5526 
-4654 
-3888 
'3224 
'2662 
-1630 
'0996 


2 *8994 
2 -2367 
1 -8982 
1 -6679 
1 -5093 
1 -3942 
1 -3083 
1 -2429 
1-1915 
1 -1434 
1 -0887 
1 -0525 



In a future paper we will show how the method can be applied to the 
problems of the motion of two spheres in a perfect fluid. 

In conclusion, I wish to express my deep sense of obligation to Dr. L» K G-. 
Filon, F.K.S., and Mr. E. Cunningham, M.A., for the valuable assistance they 
have given me, and to Prof. A. W. Porter, F.E.S., for his helpful interest and 
encouragement throughout the preparation of this paper. 



